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Abstract. In the present paper the Ising model with competing binary 
J and Ji interactions with spin values ±1, on a Cayley tree is considered. 
We study translation-invatiant Gibbs measures and corresponding free 
energies ones. 
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1. Introduction 

Nowadays the investigations of statistical mechanics on non-amenable 
graphs is a modern growing topics ([L]). One of the non-amenable graph is 
a Cayley tree. The Cayley tree is not a realistic lattice, however, its amazing 
topology makes the exact calculation of various quantities possible [L]. It is 
believed that several among its interesting thermal properties could persist 
for regular lattices, for which the exact calculation is far intractable. Here 
we mentions that in 90's a lot of research papers were devoted to studying 
of the classical the Ising model, with two spin values ±1, on such Cay- 
ley tree (see [Pr], [BG],[BRSSZ],[BRZ1],[BRZ2],[GR],[R]). In the present 
paper we extend some results of the paper [GPW], where it has been inves- 
tigated the Ising model with competing interactions, with spin values ±1, 
on a Cayley tree. Note that such models was studied extensively (see Refs. 
[MTA],[SC],[Mol],[Mo2]) since the appearance of the Vannimenus model 
(see Ref.[V]), in which the physical motivations for the urgency of study 
such models was presented. In all of these works no exact solutions of the 
phase transition problem were found, but some solutions for specific param- 
eter values were presented. On the other hand while studying such models 
it was discovered the appearance of nontrivial magnetic orderings [MR] . 

In the paper using the methods of [MR] we exactly solve a phase transition 
problem for the model, namely, we calculated critical curve such that there 
is a phase transitions above it, and a single Gibbs state is found elsewhere 
(cp.[GPW]). We also find ground states of the model. Besides, we will find 
a form of the free energy of the model under consideration. This gives us to 
study some asymptotics one. 



2. Preliminaries 

Recall that the Cayley tree of order /c > 1 is an infinite tree, i.e., a 
graph without cycles, such that each vertex of which lies on /c + 1 edges. 
Let r*^ = (y. A), where V is the set of vertices of F*^, A is the set of edges 
of r'^. The vertices x and y are called nearest neighbors, which is denoted 
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by / =< x,y > if there exists an edge connecting tlicni. A collection of 
the pairs < x,xi >,...,< X(i-i,y > is called a path from x to y. Then the 
distance d{x,y),x,y G V, on the Cayley tree, is the length of the shortest 
path from x to y. 

For the fixed x° G F we set 

Wn = {x€ V\d{x, x°) = n}, Vn = U'^^iWm, 
Ln = {I =< x,y >e L\x,y G Vn}. 

Denote = d{x,x^), x eV. 
Denote 

S{x) = {y£ Wn+i : d{x,y) = 1}, x£ Wn- 

The defined set is called the set of direct successors. Observe that any vertex 
a; 7^ x° has k direct successors and has A; + 1. 

Two vertices x,y V is called one level next-nearest-neighboring vertices 
if there is a vertex z & V such that x, y G S{z) and they are denoted by 
> x,y <. In this case the vertices x,z,y are called ternary and denoted by 
< >. 

Proposition 2.1. [G] There exists a one-to-one correspondence between 
the set V of vertices of the Cayley tree of order k > 1 and the group of 
the free products of k -\- 1 cyclic groups of the second order with generators 
ai, 02, flfc+i- 

Consider a left (resp. right) transformation shift on Gk defined as: for 
go G Gk we put 

Tg„h = goh (resp. Tp^/i = hgo,) Mh e Gk- 

It is easy to see that the set of all left (resp. right) shifts on G^ is isomorphic 
to the group G^. 

We consider models where the spin takes values in the set $ = {—1, 1} . 
A configuration o" on F is then defined as a function x G F — > (j{x) G 
the set of all configurations coincides with VL = <I>^. The Hamiltonian of the 
Ising model with competing interactions has the form 

H{a) = -J Yl ^i^Xy) - -^1 E <^)<y) (2-1) 

>x,y< <x,y> 

where J, Ji G M are coupling constants and cr G $7. 

As usually one can define Gibbs measures of the model under considera- 
tion (see, for example, [S],[MR]). 

3. On Gibbs measures 

In this section we give the construction of a special class of limiting Gibbs 
measures for the Ising model on a Cayley tree with competing interactions. 

Let /i : X — > M be a real valued function of x & V. Given n = 1,2,... 
consider the probability measure n^^^ on defined by 

/x(")(an) = Z-^ exp{-pH{an) + ^ hM^)}, (3.1) 

xeWn 
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where 

H{an) = -J ^ crn{x)an{y) - Ji ^ (rn{x)an{y), 

>x,y<:x,yeVn <x,y>:x,yeVn 

and P = ^ and (T„ : x G — ^ (T„(x) and is the corresponding partition 
function: 

The consistency condition for fi^^\an),n > 1 is 

= (3.2) 

where cr("^ = {a{x),x G Wn}- 

Let Fi C V2 C ... W^^iVn = V and //i, /X2, ■•■ be a sequence of probabihty 
measures on <I>^i, <I>^^, ... satisfying the consistency condition (3.2). Then, 
according to the Kolmogorov theorem, (see, e.g. Ref. [Sh]) there is a unique 
limit Gibbs measure /x/j on Q such that for every n = 1, 2, ... and (T„ G 
the following equality holds 

M({^|y„=^n}) =/x^")K). (3.3) 

The following statement describes conditions on hx guaranteeing the con- 
sistency condition of measures iJ,^'^\an)- In the sequel for the simplicity we 
will consider the case k = 2. 

Theorem 3.1. The measures ii'^'^\an), n = 1,2, ... satisfy the consistency 
condition (3.2) if and only if for any x eV the following equation holds: 

_i / elee^^^y+^^^ + eije^'^y + e^^-) + e \ 

here 9 = e^^*^, 9i = e'^^'^^ and < y,x, z > are ternary neighbors. 

Proof. Necessity. According to the consistency condition (3.2) we have 

Y ex.p{-PHn-li(Jn-l) + 

xeW„-i,y,z£S{x) x£W„-i,y,z€S{x) 

Y Y K^^y)^ = exp{-/3i?'„_i(cr„_i) + Yj K(^{x)}. (3.5) 

x&Wn-i yes(x) xeWn-l 
Whence we get 

n ^M(3JMx)i'T{y) + a{z))+f3Ja{y)aiz)+ 

0-(") X&Wn-l 

hy'^iy) + hz(^{z)} = Yi exp{/ia;C7(x)}. (3.6) 

xeWn-1 
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Let X G Wn-i and S{x) = {y, z}, = {cr{y), (j{z)}. Then it is easy to 

see that a^") = UxeWn-i'^x^^ • Hence 

hy(^iy) + hz(T{z)} = Yi ew{hx<y{x)}. (3.7) 

xew„-i 

Now fix X G Wn-i and rewrite (3.7) for the cases a{x) = 1 and a{x) = —1 
then we can find 

E4-)={^(j^),^(^)} exp{/3Ji(cr(y) + aiz)) + (3Ja{y)a{z) + /lyCj(y) + /i^o-(z)} 
E4")={^(2^),^(^)} exp{-/3Ji(a(y) + a(z)) + f3Ja{y)a{z) + hya{y) + /i2a(z)} 

= exp{2/t^}. (3.8) 

Denote 

W^i = exp(2Ji/3 + J/3 + + /i^) + exp(- J/3 - /ij, + /i2)+ 
exp(— J/3 + hy — hz) + cxp(— 2Ji/3 + J/3 — /i.^ — /i^) 
W-i = exp(-2Ji/3 + J/3 + /ij^ + /i^) + exp(- J/3 - /ij; + /i2)+ 
exp(— J/3 + hy — hz) + exp(2Ji/3 + Jf3 — hy — hz)- 
It then follows from (3.8) that 

exp{2/ij = (3.9) 

The equality (3.9) implies (3.4). 

Sufficiency. Now assume that (3.4) is valid, then it implies (3.9), and 
hence (3.8). Prom (3.8) we obtain 

^ exp{pji{a{y) + a{z))a{x) + pja{y)a{z) + hya{y) + hza{z)} = 

a{x) exp{a{x)hx}, (3.10) 
where a{x) = ±1. This equality implies 

Y[ J2 ex.p{pJi{a{y)+a{z))a{x)+pJa{y)a{z)+hya{y)+hza{z)} = 

^e^»-i4'''={<^(2/),a(2)} 

Yl a{x)e^p{a{x)h^}. (3.11) 

Writing A„ = rTxeWn ^i^) from (3.11) 

Zn-lAn-lfl^^-'\an-l) = Z„ //„(a„_i, a^")). 

<7(") 

Since each /i'-"-', n > 1 is a probability measure, we have 

Y: ^M(")(a„_„aW) = 1, ^ M("-^)(a„-i) = 1. 
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Therefore from these equahties we find 

Zn-lAn-l = Zn, (3-12) 

which means that (3.2) holds. This completes the proof. 

According to Theorem 3.1 the problem of describing of Gibbs measures 
is reduced to the description of solutions of the functional equation (3.4). 

By Proposition 2.1 any transformation S of the group Gk induces a shift 
automorphism 5 : f2 — > by 

{Sa){h) = a{Sh), h e Gk, a e n. 
By Qk we denote the set of all shifts of 0,. 

We say that a Gibbs measure /lx on 17 is translation - invariant if for any 
T eGk the equality n{T{A)) = fi{A) is vahd for all A G J^. 

The analysis of the solutions of (3.4) is rather tricky. It is natural to begin 
with the translation - invariant solutions where = h is constant for all 
X & V. It is clear that a Gibbs measure corresponding to this solution is 
translation-invariant. This case has been investigated in Ref.[GPW]. 

In this case from (3.4) we virtue 

where u = e^^. 

10 

Proposition 3.2 [GPW].//^i > \/3 and > then for all pairs 

{9,9i) the equation (3.13) has three positive solutions u\ < U2 < ^3, here 
u\ = l. Otherwise the eq. (3.13) has a unique solution n* = 1. 

Remark. The numbers u\ and u\ are the solutions of the following 
equation 

ii2 + (l + a)M + l = 0, (3.14) 

here a = — 6\. Hence u\u'^ = 1 and if /3 ^ cx) then U3 ^ 00 and u\ 0. 

By /Lti, /X2, Hz we denote Gibbs measures corresponding to these solutions. 
Denote Ux = exp(2/i3;),x G V. Then the functional equation (3.4) is 
rewritten as follows 

"^"^ 6uyu, + Oiiuy + u,) + ele ^ ' ' 

here as before < y, x, z > arc ternary neighboring vertices. 

201 

Proposition 3.3. Let 61 > \/3, 9 > — - and Ux be a solution of 

9i — 3 

equation (3.15). Then 

Ui ^ Ux ^ for any x & V. 



The Proof is similar to the proof of Proposition 4.3. of [MR]. 
From this proposition we infer the following (See [Ge]) 
Theorem 3.4. For the model (2.1) with parameters Ji > and J G M 
on the Cayley tree the following assertions hold 
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(i) if 6*1 > \/3, 6 > — - then the measures jii and ^3 are extreme; 

(ii) in the opposite case there is a Gibbs measure n*{= ^12) and it is ex- 
treme. 

Remark. This theorem specifies the result obtained in [GPW] as they 
has proved that a phase transition occurs if and only if the above indicated 
conditions is satisfied, and the extremity was open. The formulated theorem 
answers that the found Gibbs measures are extreme. Note that in the (i) 
case the measure 1^2 need not to be extreme, some relevant information 
concerning the extremity of 112 for the Ising model can be found in [BRZl]. 

It is clear from the construction of the Gibbs measures that the measures 
jii and Hz depend on parameter (3. Now we are interested on the beheviour 
of these measures when (3 goes to 00. 

Put 

(T_|_ = {(t{x) : (t{x) = 1, Vx € r^}, 

(7_ = {a{x) : a{x) = -1, Vx G T^}. 

29 

Theorem 3.5. Let 6\ > \/3 and 6 > , then 

el- 3 

fJ-l Sa_ , ^3 as /? ^ 00, 

here 5^ is a delta-measure concentrated on a. 

Proof. Consider the measure /X3. This measure corresponds to the func- 
tion hx = h3, X e V, here ^3 > (see Proposition 3.2). Let us first consider 
a case: 

fisiaix) = 1) = , = -^j— ^1 as /3 ^ 00, 

since U3 ^ 00 as /3 ^ 00, here x £ V. Let us turn to the general case. From 
the condition imposed in the Theorem we find that Ji > 0. Now separately 
consider two cases. 

First case. Let J > 0. Then from the form of Hamiltonian (2.1) it is 
easy to check that H{an\vn) ^ ^i'^+lvn) for all a G $7 and ra > 0. It follows 
that 

„ . I ^ eM-PHia+\vJ + h3\Wn\} 

M3(0-+|\/„) ^ 



E exp{-/3i?(a„) + /i3 E ^(x)} 

1 

exp{-PH{an) + h3 Yl ^i^)} 

exp{-PHia+\vJ + h3\Wn\} 



> 



> r, — ^ , ^ — 1 as /3 — 00. 



The last inequality yields that /X3 — (5, 



0-+- 
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Second case . Let J < 0. Let us introduce some notations. 
>x,y<:x,yeVn 

<x,y>:x,y£V„ 

CK)= ^a(x), C = C{a+\vJ. 

xeVn 

Then it is easy to see that the following equality holds 

M3(C^+|V„) = 1 

We want to show that 



(,jp{A-A{an))eJiP{B-B{an))eh3{C-C{ar,)) 



X] gJ/3(A-A(a„))gJi/3(B-B(a„))g?i3{C-C(a„)) ^ ^ ^ P ^ OO. 

It is enough to prove that 

gJ/3(A-A(a„))gJi/3(B-B(CTn))gft3(C-C(CTn)) ^ ^ ^ ^ OO 

for all an € fiy„, cr„ 7^ (7+|y„. We rewrite the last sentence as follows 
1 1 



gJ/3(A-A(S„))gJi/3(£;-B(S„))gft3(C-C(S„)) g,(A-A(5„))/2^p--B('^n))/2^^^^(c_c((T„))/S 

(^2 _ 3)(A-A(ct„))/2 

^ ^(A-A(5„)+B-B(a„))/2^^,^(C_c(a„))/2 ~ 
(02 _ o\iA-Aia„))/2 
- ^(A-A(a-„)+B-B(5-„))/2 . • V^--""^ 



< 



'1 



here we have used the inequality > If^ . 
Obviously, if /3 is large enough we have 

(5/2 _ 3^(A-A(5„)) g|2(A-A(a„)) 



QA-A{a„)+B-B{an) QA-A{an)+B-B(a„) 



IQ2 _ 3)(A-A(5„)) 

B{an)—A{an) < B—A then the last relation implies that 



jA-A{d-„)+B-B{a„) 

is bounded, and hence from (3.16) we get the desired relation. 
Now it remains to prove the following 

Lemma 3.6. For every n > and an G Jly„ the following inequality 
holds 

B{an) - A{an) < B - A. (3.17) 
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Proof. Denote C(ct„) = {x & Vn ■ = — !}• Maximal conneeted 

components of C((7„) we will denote by /Ci(<Jn)) ' ' ' ) ^m(c"n)- For a connected 
subset K. oiVn put 

dK = {x e Vn\]C : < x,y > for some y e /C}. 

= {x G y„ \ /C : > X, y < for some y G /C}. 
Prom definition of A{an) and B{an) we get 

i?(a„) = i?-2^|a/C,K)|, 

A{an) = A - 2^ \ U^^jlCm{an)\, 

j 

here |A| stands for a number of elements of a set A. 

To prove (3.17) it enough to show that |3^/C| < |9/C| for all connected 
subsets /C of Vn- For each a; G we can show a y = y{x) G Indeed, 
a > x,t <, t E K. and < a;, > then y{x) = x if z G /C and ?/(x) = z if 
z ^ fC. It is clear y(x) G dIC. Now we will prove that if xi 7^ X2 G d^fC then 
y{xi) 7^ y(x2) G 5/C. Let < xi,zi,ti >,< X2, 2:2,^2 >, where ti,t2 G /C. By 
definition of y{x) we have y{xi) G {xj, z,}, z = 1, 2. So to prove y(xi) 7^ y{x2) 
for xi 7^ X2 it is enough to show zi ^ Z2- Note that in our case (i.e. k = 2) 
if xi 7^ X2 G d^KL then d{xi^X2) > 3, since < Xi,Zi >,i = 1,2 and hence we 
get ^1 7^ 2:2- Thus |3^/C| < \d)C\. So Lemma is proved. 

By similar argument the theorem can be proved for the measure fii . Thus 
the theorem is proved. 

From Theorem 3.5 we conclude that cr+ and (T_ are ground states of the 
considered model. 

Remark. If in the condition of Theorem 3.5, we put J = then the ob- 
tained result coincides with Theorem 2.3 of [BRZ2]. From proved theorem 
we infer that the measure 112 is not extreme while /3 is large. For small f3 
the extremity of fi2 would be considered elsewhere. 



4. A FORMULA OF THE FREE ENERGY 

Consider the partition function Zn{l3,h) of the state /x^ (which corre- 
sponds to solution h = {hx,x G V} of the equation (3.4)) 

Zn{P,h)= J2 exp{-/3i7(a„) + ^ /i^a(x)}. 
The free energy is defined as 

The goal of this section is to prove following 
Theorem 4.1. 
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(i) The free energy exists for all h, and is given by the formula 

F(/3,^) = - lim ^lim ^ Dp{J,, J,hy,h,), (4.2) 

where y = y{x),z = z{x) are direct successors of x; 
D{Ji, J, hy, h^) = djp{Ji(i + h,) + djp{-,h[3 + h^) + 
^jAh + f{-Ji(i + K; 6)- hy + f{.hP + h,- 9))- 

dp{x) = ln[4cosh(a; — (3) cosh(x + /?)]; 

^pix^y) = ^ln[4cosh(x-/3)cosh(y + /3)]; (4.3) 

f{x,e) = tanh"^(6'tanhx), 9 = tanh(J/3). 

(ii) For any solution h = {hx^x G V} of (3.4) 

F{(5,h) = F{f5,-h), 
where —h = {—hx, x G V}. 
Proof, (i) We shall use the recursive equation (3.12): 

— -^n—l^n—li 

where An = Y\ ^(^) a{x) = a{x, Ji, J, P),x G F is some function, 

x£W„ 

which we will define below. Using (3.10) we have 

a{x) = 4 v/cosh(Ji/3 - J/3 + h^) cosh(Ji/3 + + h^) x 
x/cosh(- Ji^ - J/3 + h^) cosh(- Ji/3 + J/3 + h^) x 

^cosh(- Ji/3 + hy + f{-JiP + h,; 9)) cosh(Ji/3 + hy + /(Ji/3 + /i^; 9)). 
Thus, the recursive equation (3.12) has the following form 

ZniP;h)=exp( ^ Df3{Ji,J,hy,h,)]Zn-iiP,h). (4.4) 

This gives (4.2). Now we prove existence of the RHS limit of (4.2). By 
proposition 3.4 we have hx G Inu^; ^ lnu3] consequently functions from 
(4.3) and so D is bounded i.e. |D(Ji, J, /iy, ^z)| < for all hy,hz- Hence 
we get 



1 

7^ E E Dp{J,,J,hy,h,)< 



3-2 

k=l+l X&Wn-k 

^ 2-'=-i < Cfs ■ 2-K (4.5) 

fc=Z+l 

Therefore, from (4.5) we get the existence of the limit at RHS of (4.2). 
(ii) Now, we shall prove that 

F{P,-h)=F{P,h). (4.6) 
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It is easy to see that if h = {hx,x G V} is a solution to (3.4) then —h = 
{—hx,x G V} also is solution to (3.4). The equality (4.6) follows from the 
following equality 

D{Ji, J, -hy, -hz) = D{Ji, J, hy, hz), 

which is a consequence of the following properties 

1) dp{-x) = dp{x); 

2) f{-x-e) = -f{x-e)- 

3) /:^p{-x,-y) = /\p{y,x). 
The theorem is proved. 

The rest of the section is devoted to study of asymptotical properties of 
the free energy F(/3, /i) as /? ^ oo for any constant h i.e. hx = const, x &V. 
In this setting F(/3, h) has the form: 

F{P,h)=D{Ji,J,h,h). (4.7) 

Denote 

F(oo) = lim F{P,h). (4.8) 

By Proposition 3.2, we know that the equation (3.4) has exactly three 
constant (translation -invariant) solutions: h = h{(3) = {hx = ^lnul,x G 
V}, {hx = 0,x e V} and {hx = ^lnul,x G V}. Using (3.14) one can find 
that h has the following asymptotic 

h{P) = M(3 + o{p-^), N>2, as P^oo (4.9) 

where M = ^ max{2(Ji - J); 3Ji - J; 4Ji; Ji - J, 0}. Using (4.9) it is easy 
to see that 

djpiJip ±h) = ^Y.\Ji+eJ±M\ + o{p-^), 

e=±l 

f{±Ji(3 + Mp + o{(3-''y,e) = ^(|±Ji- J + M|-|±Ji + J + M|) + o(/3-^); 
here we have used the following easily checking asymptotic 

In (2 cosh(a/3 + o(/3-^))^ = (3{\a\ + o{(3-^)), 

where N >2 and a G M. Consequently, 

F(/3,/i)=F(oo) + o(/3-^), (4.10) 

where 

F(oo)= ^ 5Ji+M-^ Yl \SJi+PJ + M\e + \Y1 \Ji+£J + SM\] . 
s=±i \ £=±1 £=±1 / 
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5. Discussion of results 

It is known that to exact calculations in statistical mechanics are paid 
attention by many of researchers, because those are important not only for 
their own interest but also for some deeper understanding of the critical 
properties of spin systems which are not obtained form approximations. So, 
those are very useful for testing the credibility and efficiency of any new 
method or approximation before it is applied to more complicated spin sys- 
tems. In this paper we have exactly solved an Ising model on a Cayley tree, 
the Hamiltonian of which contains the nearest-neighbor and competing in- 
teractions, namely, we calculated critical curve such that there is a phase 
transitions above it, and a single Gibbs state is found elsewhere. It is found 
ground states of the model. Wc also explicitly express the free energy of the 
model under consideration, associated with the translation invariant Gibbs 
measures, these enable us to find some asymptotics of the free energy as 
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